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Abstract. The present paper is devoted to find a new generalization of the generalized 
Chaplygin gas. Therefore, starting from the Hubble parameter associated to the Chaplygin 
scalar field and using some elliptic identities, the elliptic generalization is straightforward. 
Thus, all relevant quantities that drive inflation are calculated exactly. Finally, using the 
measurement on inflation from the Planck 2015 results, observational constraints on the 
parameters are given. 
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1 Introduction 

Clearly, evidence from the most varied observations suggest that near 25% of the cosmic 
energy-matter corresponds to the so-called cold dark matter (CDM) and the remaining 75% 
corresponds to the so-called dark energy (DE). The dark matter problem born with the 
Zwicky’s observation of a large velocity dispersion of the components in the Coma cluster [1] 
together with Babcock observation of galactic rotation curves in the Andromeda galaxy [2]. 
This paradigm has been studied in several context, for example from the particle physics point 
of view, the candidates are axions, inert Higgs doublet, sterile neutrinos, super-symmetric 
particles and Kaluza-Klein particles [3-7], while from the geometrical point of view Weyl’s 
theory has received much attention from theoretical and experimental physicists [8-15]. On 
the other hand, one of the most natural candidates to explain DE corresponds to a cosmo¬ 
logical constant A, which, although successful as a fit model, it shows a huge discrepancy 
between its theoretical and observed values [16]. This evidence has promoted the study of 
new types of candidates that can link the different stage of the universe. One type of cosmic 
fluid is often called quintessence and has been widely used because it behaves like a cosmo¬ 
logical constant by combining positive energy density and negative pressure. Nevertheless, 
unlike the cosmological constant, the pressure of this fluid is dynamic. An example of this 
type is the so-called Chaplygin gas (CG) which is an exotic fluid having the equation of state 
(EoS) 

B , , 

Pch = -, (1-1) 

Pch 

where B is a positive constant, pch and pch are the pressure and energy density in a comoving 
reference frame, respectively. This EoS was originally proposed to describe the lifting force 
on a wing of an airplane in aerodynamics [17]. In cosmology, it was first introduce by [18-20], 
nevertheless, in its original form, this fluid presents inconsistency with some observational 
data [21-23] such that a generalization is necessary to improve consistency with data. A first 
attempt was made considering the generalized Chaplygin gas (GCG) whose EoS is given by 
[24-27] 

B , , 

Pgcg — “ 77 ^’ 

Pgcg 

where B > 0 and the GGG parameter lies in the range 0 < a < 1. Obviously, the GG is 
recovered by making a = 1, while the case a = 0 mimics the effect of a cosmological constant. 
An important characteristic of the Ghaplygin gas (GG and GGG) is that this interpolates 
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between a dust dominated phase where p oc a~^, and a de Sitter phase where p = —p. 
The different parameters that enter into the model has been confronted successfully with 
many observational data (see, for example [28-30]), but, unfortunately, the model produce 
oscillations or exponential blowup in the matter power spectrum, being inconsistent with 
observation [31]. 

On the other hand, an approach to single-field inflation using the GCG as inflaton 
field was performed by del Campo [32]. By comparing with the measurement for the scalar 
spectral index together with its running from the Planck 2013 data [34], del Gampo found 
that the best value for the a-parameter is given by a = 0.2578 ±0.0009. Therefore, following 
the idea outlined in [35], the main goal of this paper is to give a new generalization for the 
Chaplygin scalar field by using elliptic functions to describe the inflationary stage and then 
perform the confrontation between the model and the measurement recently released by the 
Planck 2015 data. 

The outline of this work is as follows: the basic aspects to describe exact solutions 
to inflationary universe models within the framework of the Hamilton-Jacobi approach to 
cosmology, together with review of the results found earlier by del Campo for the GCG are 
showed in Section II. In Section III the Chaplygin-Jacobi gas is presented and the Hamilton- 
Jacobi formalism is applied to find the most relevant quantities in the inflationary stage. 
Finally, in Section IV conclusions and final remarks are presented. In addition, Appendix A 
few useful formulas and identities that satisfy the Jacobi elliptic functions are shown. 


2 The generalized Chaplygin gas in inflation 


The fundamental quantity in inflation is the scalar inflaton field (p which leads to the universe 
to expand extremely rapid in a very short time. The evolution of this field becomes governed 
by its scalar potential, V{(j)), via the Klein-Gordon equation 

(f) ± 3H(j) ± = 0, (2-1) 


where represents a derivative with respect to cp. Thus, this equation of motion, together 
with the Friedmann equation. 




Stt 

3m| 


\^‘' + V{cP) 


( 2 . 2 ) 


obtained from Einstein general relativity theory, form the most simple set of field equations, 
which could be applied to obtain inflationary solutions. The set of equations (2.1-2. 2) are 
the mainstay of Hamilton-Jacobi approach [36-42]. 

In order to obtain exact solutions in this approach, del Campo [32] and Dinda et al [33] 
used the Chaplygin gas whose generating function is given by 


H{(p) = Hq coshi+“ [(1 ± a)<I>)] , 


(2.3) 


where 4? = ^/^{(p — (po) is a dimensionless scalar field and Hq, the value of H 

is a constant given by Hq = H2(i+“). From the generating function 

obtain all the relevant quantities for inflation, for example the scale factor 


when (p = (pQ, 
(2.3) one can 


a{(p) = Oi exp 


47r H{(p) 

mj J^. H^^{(p) 
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(2.4) 










which becomes ^ 




( sinh[(l + a)<h] \ 
Vsinh[(l + a)^i]) 


2 

3(l+a) 


(2.5) 


The number of e-folding of physical expansion that occur in the inflationary stage is given 

by 



dvr H 

2 / tr ) 

J (pf 


( 2 . 6 ) 


so, del Campo found that for the GCG this quantity becomes^ 


iV(0) = 


sinh[(l + a)<he] 
3(1 + a) \ sinh[(l -t- 0 ;)$] 


■ In 


while the scalar potential 

V{cj>) = 

is given by the following expression: 


/ 3ml\ 



^- -Hi 

Utt 



cosh i+“ [(l + a)<h] 


(2.7) 


( 2 . 8 ) 


(2.9) 


-t j- 

where Vq = . In the slow-rolls approximation this scalar potential is well behaved and 

reduces to 

14 _r (<(>) —Vo cosh 1 +^ [(1 + a)<h]. (2.10) 

In next section, the Hamilton-Jacobi formalism is applied to obtain exact solutions in 
the inflationary stage using a elliptic generalization of the GGG. 


3 The generalized Chaplygin—Jacobi gas 

It is well known that trigonometric and hyperbolic functions are special cases of more general 
functions, the so-called elliptic functions, which are doubly periodic in the complex plane, 
while the above are periodic on the real axis and the imaginary axis, respectively. With this 
in mind and using the rules given in appendix A, the generating function (2.3) is conveniently 
written as 

H{(j),k) = Ho nc^ [{1 +a) (3.1) 

where nc(x) = l/cn(x), and cn(x) = cn{x\k) is the Jacobi elliptic cosine function, and 
k is the modulus. This choice allows to obtain the generating function of the GCG (2.3) 
making A: —?■ 1 in equation (3.1). Henceforth we denote this scalar field as the generalized 
Chaplygin-Jacobi gas (GCJG). 

Taking into account that pressure and energy density are given by 

(3. 


'^the original paper [32] presents a typo in this formulae, 
^the original paper [32] presents a typo in this formulae. 
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Figure 1: Plot of the pressure against density for generalized Chaplygin-Jacobi gas (3.6). 
Here we plot different values of the modulus k and GCG parameter a. In all graphs was 
used B = 1. 


and 

3m^ 9 

P<p = (3.3) 

respectively, and by using the generating function (3.1) it is no hard to obtain that 


_B2(l+a) 2(2+a) 

Pci> = --nc i+“ [(1 + a) (dn^ [(1 + a) $] - A:), (3.4) 

and ^ 

/90 = i?2(i+“) ncH^ [(1 + a) $] , (3.5) 

where k' = 1 — k \s the complementary modulus. Therefore, by using eq. (3.5) one obtain 
nc [(1 + a) <1>] as a function of and then, by substituting into eq. (3.4) together with some 
identities of the Jacobi elliptic functions, the EoS related to the GCJG read 

= + (3.6) 


Obviously, in the limit A: —)■ 1 (or k' —)■ 0) this EoS coincide with ( 1 . 2 ). Note that, for A: 7 ^ 1, 
the GGJG presents a novel behaviour in which the pressure is a positive quantity. This fact 

(c) (c) 

occurs at the range py < p^, where critical density, , depends on the parameters across 
the following expression 



B 



1 

1+a 


(3.7) 
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At this point = 0, and the fluid presents a dust-like behaviour. An interesting feature of 
this EoS appears when A: = 0 (or k' = 1) which corresponds to the trigonometric limit. In 

addition to the root obtained from equation (3.7), must also add = 0, and therefore the 

1 

GCJG presents a dust-like behaviour when Pcj, = 0 and p^ = (2 B) i+“ , as shown in Figure 1. 
For this fluid one has simply that the EoS parameter becomes 


= 


P± 

P<f> 


^ l+a 


^-2k' 

1+a 

P<t> 


(3.8) 


while the square of the speed of sound is 


_2 dp^ (2 + a)/c' ^, 


(3.9) 


Note that the speed of sound vanished when the density reaches the value 



B 



— pQ^(2 + a) ^ 

1 

2 + q; I 


1 

1 + a 


(3.10) 


Of course, the existence of these points is restricted to pairs that satisfy the condition k~^ > 
(1 + a)^. For example, if a = 0, then k < 1, while for the case a = 1 the modulus must be 
k < 0.25. It is important to note that null speed of sound do not implies null pressure. In 
fact, at this point the pressure is 

Ppo = P^ipp = OK + ^(1 + «) ± \/^' [1 - k{l + a)^]\. (3.11) 

(2 + a)2(4 )“ 

On the other hand, using eq. (2.8) together with eq. (3.1) one obtain for the scalar 
potential 

^ 2 nc^ [(1 + g) $] - dc^ [(1 + a) $] sc^ [(1 + a) $] 

nci+“ [(1 + a) $] 

1 1 

where Vq = ii? 2 (i+a). This scalar potential is depicted in Figure 2 for different values of the 
GCG parameter a and the modulus k. Also, by using eqs. (2.8) and (3.1) is straightforward 
to found the potential as a function of the generating function, so 


V{H) 



k' 


U) 




2(l+a) 


(3.13) 


By combining the generating function (3.1) with eq. (2.4) allows to find the scale factor 
as a function of the inflaton field: 


a{(t)) = Qi 


sd[(l + a)<h] \ 3(1+0) 

sd[(l + a)4>i]y 


(3.14) 


Here sd(x) = sd(x| A:) = sn(x)/dn(x), where sn(x) = sn(x| A:) is the Jacobi elliptic sine 
function and dn(x) = dn(x| k) is the Jacobi elliptic delta function. 
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Figure 2: Plots of the scalar potentials for generalized Chaplygin-Jacobi gas, F($), as a 
function of the dimensionless scalar field, $ = ^^{(f) — <^o)- Here the scalar potential V (‘h) 


is expressed as a multiple of the constant Vq = 


With the aim to verify that an inflationary period occurs, it is very helpful to introduce 
the deceleration parameter Q = -^ [ ] — 1, which becomes 


3 dn^ [(1 + a) <I>] sn^ [(1 + a) d>] 

^^2 cn 2 [(1 + a ) $] ■ 


(3.15) 


Making use of some elliptic identities it is possible obtain the value for (jig, the value for the 
scalar held when inflation ends (i. e. g(i^e) = 0). For example, if eq. (3.15) is writes in terms 
of the Jacobi elliptic sine one arrives to the following quadratic equation 


3fey2-5y + 2 = 0, 


(3.16) 


where y = sn^ [(1 + a) <l>e]. Therefore, solving for <he one get that the four roots are given by 


<f>e = 


F [arcsin {F^/u±) , k] 


1 “h Oi 


(3.17) 
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where F{ip, k) is the normal elliptic integral of the hrst kind and 



(3.18) 


It is important to note that real solutions are found for all 0 < A: < 1. Also, the condition 
y < 1 implies that the allowed solution are -iiyfyZ. 

On the other hand, using eq. (2.6) one get that the number of the e-folds becomes 


1 V( 0 ) 


2 f sd[(l + «)che] l 
3(1 +a) 1 sd[(l + a)$] /’ 


where <he is given by eq. (3.17). For example, for = 70 the above leads to 


(3.19) 


<I>70 — 


arcsm 


■s/l+fca 


1 “1“ Q! 


(3.20) 


where x is given by 

X = sd[(l + a)^e] . (3.21) 

In order to obtain some information about the perturbations, result instructive to dehne 
the three first hierarchy parameters in terms of the generating function and its derivatives. 
Thus, the first Hubble hierarchy parameter en is given by 


dhiH 

din a (dvryyFI/ 


(3.22) 


Note that this parameter provides information about the acceleration of the universe, so 
during inflation the bound e// <C 1 is fulhlled, and inflation ends once e// = 1. Taking 
account that en = q + 1-, one obtain from (3.15) that 


_ 3 dn^ [(1 + a) <I>] sn^ [(1 + a) <I>] 
^ 2 cn 2 [(1 + a ) ^>1 

The second Hubble hierarchy rjH is defined as 

~ dlna (dvr J V ^ y ’ 


(3.23) 


(3.24) 


where plugging the Chaplygin-Jacobi scalar field becomes 


1 + a cn2[(l + a)<h] {l + a)k' 


+ 


sn2[(l + q:) 4)] ' dn^[(l + a)<I>] 

Finally, the third Hubble hierarchy parameter, is dehned by 


(3.25) 






(j) ^ ,4>4^4> 

ip 


-7 ~ 


(3.26) 


















which becomes 




(2a^ + 7a + 6) k' 


sn^ [(1 + a)$] dm [(1 + a)<h] 


+ [{2k - 1) - ka{l + 2a) cn^ [(1 + a)<h]] | 

J 


(3.27) 

The quantum fluctuations produce a power spectrum of scalar density fluctuations of the 
form [43, 44] 


ff 


rn(k) = l-l (-1 


ff\- 


(3.28) 


aH=k 


This perturbation is evaluated when a given mode k crosses outside the horizon during 
inflation, i. e. at aH = k. These modes do not evolve outside the horizon, so we can assume 
they keep fixed value after crossing the horizon during inflation. In order to obtain some 
comparison with the available observational data, we introduce the scalar spectral index Ug 
dehned as 

dlnT^T^ 


n,o - 1 = 


dink 


After a brief calculation, one obtains that 

ris - 1 = -4e/f + 2r]H, 

so, using Eqs. (3.23) and (3.25), the spectral index becomes 

a k' 


n, — 1 = 3 


— {2k — 1) + {2 + a)k cm [(1 + a)4>] > . 


(3.29) 


(3.30) 


(3.31) 


cn^ [(1 4- a)$] 

As a first observation, note that the above equation allows to write cn^ [(l+a)^^] as a function 
of Ug. By solving the quadratic equation on cn^ [(1 + a)^^] is straightforward to find that 

i(n, -l) + {2k-l) \ 4A:A:'a(2 + a) 'I 


C{ng, a, k) = cm [(1 + a)^] = 


2 k {2 +a) 


[\{ng - 1) + {2k - l)]^ ] ■ 
(3.32) 

From here, it is no hard to see that Ug get the value equal to one when the dimensionless 
scalar field get the value 

F [arcsin (±Va±) ’ 


where 


X± = 


4*718 = 1 ~ 


2k - 1 
2k{2 + a) 


1 T a 


i± Wi- 


4:kk' a {2 + a) 


( 2 /c- 1)2 


(3.33) 


(3.34) 


Ur 


On the other hand, one can introduce the running scalar spectral index defined by 
_ which results to be 


dink 


Ur 


— 10 ch rjH — Se^jj — 2 . 


(3.35) 


Therefore, using Eqs. (3.23), (3.25) and (3.27) one get that 

= —6(1 + a) {fc' a nc^ [(1 + a)<h] — A: (2 + a) cn^ [(1 + a)<i>]} en- (3.36) 


Ur 
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Figure 3: Plot of the transcendental relation between the a parameter and the modulus 
k for the generalized Chaplygin-Jacobi gas (3.37). The shaded region corresponds to pairs 
that are permitted by the observational data, including the value found by del Campo using 
k = 1, a = 0.2587 ± 0.0009 [32]. In this graph has been used nmn = —0.003 ± 0.007 and 
Us = 0.968 ± 0.006 from the Planck collaboration [45]. 


Note that using eq. (3.32) into eq. (3.36) one can obtain the running scalar spectral index 
as a function of scalar spectral index 


^run — 


9(1 + a)(l — C){k' + kC)[k' a — k{2 + a)C^] 


(3.37) 


and then, giving their value from the Planck collaboration [45], Umn = —0.003 ± 0.007 and 
Us = 0.968 ± 0.006, it is possible to obtain a region of validity for the transcendental relation 
between a and the modulus k which is showed in FIG. 3. Based in this results, the cases 
A: ~ 0 and A; ~ 1 are favoured in all range 0 < a < 1. The situation is completely different 
when the modulus k tends to 0.5 where small values of the GCG parameter are favoured. 

It is known that not only scalar curvature perturbations are generated during infla¬ 
tion. In addition quantum fluctuations generate transverse-traceless tensor perturbations 
[46], which do not couple to matter. Therefore they are only determined by the dynamics 
of the background metric. The two independent polarizations evolve like minimally coupled 
massless fields with spectrum 


Vt 


Ibvr 



(3.38) 


In the same way as the scalar perturbations, it is possible to introduce the gravitational wave 
spectral index ht defined by ut = becomes ut = — 2e// . At this point, it is 
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Figure 4: Plot of the parameter r as a function of the scalar spectral index Hs for the 
generalized Chaplygin-Jacobi gas (3.41). The marginalized joint 68% and 95% confidence 
level regions using Planck TT + low P, Planck TT, TE, EE + low P and Planck 2013 data 
release with k = 0.05/i Mpc“^ are shown [45]. To perform a consistent comparison, the values 
of r have been normalized in such a way that r = 0 when n* = 1, and r = 16e// is used. Note 
that this is four times the normalization used in this work, see equation (3.39). 


possible to introduce the tensor-to-scalar amplitude ratio r = ^ 


which becomes 


r = Aen, 


(3.39) 


therefore one arrives to 

dn^ [(1 + a) 4>] sn^ [(1 + a) 4>] 

cn^ [(1 + a) <!)] 


(3.40) 


With the help of eq. (3.32) is possible to find the relation between the tensor-to-scalar 
amplitude ratio and the scalar spectral index, which reads 


r{ns, 


a,k) = 6 


{k' + kC){i-C) 

c 


(3.41) 


which is contracted with the measurement from Planck collaboration [45] as is showed in 
Figure 4. 


4 Final remarks 

This paper presents an inflationary universe model in which the inflaton field is characterized 
by an EoS corresponding to a generalized Chaplygin-Jacobi gas (GCJG), i.e. — 

2 k' ^ , where a is the GCG parameter, and was considered to lie in the range 

0 < a < 1, and the k parameter is the modulus of the elliptic function which lie in the 
range 0 < A; < 1. In this study the kinematical evolution was described using the Hubble 
parameter given by H{(j),k) = [(1 + «)<!>]. From here, all relevant quantities that 
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describe inflation (scale factor, number of e-folds, deceleration parameter, etc.) were obtained 
in terms of Jacobian elliptic functions, allowing express them in an analytical form. 

As a direct application of these analytical results, It is possible to obtain an expression 
relating the running scalar spectral index with the scalar spectral index (cf. equation (3.37)). 
Using PZancfc 2015 data [45], nmn = —0.003±0.007 and Ug = 0.968±0.006, the transcendental 
relation in the a-k space, showed in Figure 4, was obtained. Notice del Campo [32] obtained 
just one value a = 0.2587 ± 0.0009 using the Planck 2013 data [34], which is in agreement 
with the present results. In the same way, a relation between the tensor-to-scalar ratio 
parameter, r, and Ug is found (cf. equation (3.41)), and contrasted with the marginalized 
joint 68% and 95% conhdence level regions using Planck TTj-low P, Planck TT, TE, EE -|- 
low P and Planck 2013 release with k = 0.05/i Mpc“^ [45]. Based on these comparisons, it 
can be argued that the GCJG model is very auspicious regarding the available observational 
tests. 

Finally, since this generalization introduces a new parameter (the modulus k), it provides 
a viable opportunity to improve confrontations with the observational data. Thus, this 
model offers new opportunities to further investigate on the accelerated stage and interacting 
models. 


A A brief review of Jacobian elliptic functions 


As starting point, let us consider the elliptic integral [47-51] 

ry dt r dO 


u{y,k) = u = 


'o Y^(l - t2)(i _ kt'^) 


f 


y/l — k sin^ 9 


= F{ip,k), 


(A.l) 


where F{(p,k) is the normal elliptic integral of the first kind, and k is the modulus. The 
problem of the inversion of this integral was studied and solved by Abel and Jacobi, and 
leads to the inverse function defined hy y = sin(^ = sn{u,k) with ip = amu, and are called 
Jacobi elliptic sine u and amplitude u. 

The function sn u is an odd elliptic function of order two. It possesses a simple pole 
of residue 1/k at every point congruent to iK' (mod 4A, 2iK') and a simple pole of residue 
— 1/k at points congruent to 2K + iK’ (mod 4 K, 2iK'), where K = K{k) = F{Tr/2,k) 
is the complete elliptic integral of the first kind, K' = F{'k/ 2,k'), and k' = 1 — /c is the 
complementary modulus. 

Two other functions can then dehned by cn(u, k) = sjl — y"^ = cos p, which is called the 
Jacobi elliptic cosine u, and is an even function of order two; dn(n, k) = -s/l — ky'^ = /S.p = 
sin p, called the Jacobi elliptic delta u which is an even function. The set of functions 
{snu, cnu, dnu} are called Jacobian elliptic functions, and take the following special values 


sn(u|0) = sintt, sn(tt|l) = tanhu, (A.2) 

cn(u|0) = cosu, cn(tt|l) = sechu, (A.3) 

dn(tt|0) = 1, dn(u|l) = sechu. (A.4) 


The quotients and reciprocal of {snu, cnu, dnu} are designated in Glaisher’s notation by 


nsu = 

1 

cs u = 

cnu 

ds u = 

dnu 

(A.5) 

snu’ 

snu ’ 

snu ’ 

ncu = 

1 

sc u = 

snu 

dcu = 

dnu 

(A.6) 

cnu’ 

cnu’ 

cnu ’ 

ndu = 

1 

sdu = 

snu 

cd u = 

cnu 

(A.7) 

dnu’ 

dnu’ 

dnu 
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Therefore, in all, we have twelve Jacobian elliptic functions. Finally, some useful fundamental 
relations between Jacobian elliptic functions are 


o 9 

sn ri J- cn u = 1, 

(A.8) 

dT?u J- k sn^u = 1, 

(A.9) 

dn^u — k cn^u = k', 

(A.IO) 

2 1 7 ^ 2 J 2 

cn u k sn u = an u. 

(A.ll) 
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